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土 中 水 膜 厚 度 变化 规律 及 未 冻 水 含量 预测 方法 


万 旭 升 ， RE =, 路 建国 ， E Mf 
(西南 石油 大 学 土木 工程 与 测绘 学 院 , 四 川 成 都 610500) 


摘 要 : 未 冻 水 含量 改变 直接 影响 冻 土 的 热力 学 及 变形 特性 ,也 是 水 热 数值 计算 的 关键 条 件 。 为 了 探究 未 冻 水 含 


量变 化 规律 ,首先 采用 等 效 粒 径 表 征 土 颗粒 粒 径 分 布 ，; 
堆积 和 立方 最 密 堆 积 平均 的 堆积 方式 ,基于 多 孔 介质 预 


a 


各 士 体 简 化 成 等 效 粒 径 球 状 堆积 体系 。 其 次 考虑 简单 立方 
融 理 论 ,提出 土 中 未 冻 水 含量 的 计算 方法 ,并 通过 试验 数据 


验证 其 合理 性 。 最 后 分 析 了 杂质 密度 以 及 有 效 粒 径 大 小 对 水 膜 厚度 以 及 未 冻 水 含量 的 影响 。 研 究 结果 表明 : 土 颗 


粒 表面 电荷 密度 引发 的 自由 能 对 水 膜 厚度 的 改变 极 易 受 杂 质 浓度 的 影响 , 随 着 浓度 增 大 , 土 颗粒 表面 电荷 对 水 膜 
的 影响 越 来 越 小 。 土 中 液态 水 含量 主要 由 土 颗粒 表面 水 膜 厚度 变化 来 决定 , 当 等 效 颗 粒 减 小 时 ,通过 间隙 水 求解 
的 未 冻 水 体积 比例 逐渐 增 大 ,并 且 在 较 低 浓度 时 发 挥 更 大 作用 。 对 于 较 大 颗粒 土 ,未 冻 水 含量 模型 预测 效果 好 。 
当 粉 质 黏 士 黄土. 砂 士 等 效 粒 径 分 别 缩小 0.28 .0.3 .0.36 倍 时 ,未 冻 水 体积 含量 计算 结果 与 试验 值 吻合 较 好 。 

关键 词 : 预 融 ; 未 冻 水 含量 ; 等 效 粒 径 ; 杂质 浓度 ; 表面 电荷 密度 


我 国 冻 土 面积 分 布 广泛 ,多 年 冻 土 和 季节 性 冻 
土 面 积 之 和 超过 国土 面积 的 70%"。 寒 区 工程 建设 
面临 因 水 热 相互 作用 引发 的 冻 胀 融 沉 现象 ,严重 制 
约 着 公路 .铁路 以 及 构筑 物 等 基础 设施 建设 。 水 分 
是 土 体 发 生 冻 融 的 物质 基础 ,未 冻 水 含量 的 改变 直 
接 影 响 冻 土 的 热力 学 及 变形 特性 。 同 时 ,未 冻 水 含 
量 随 温度 变化 规律 也 是 水 热 数值 计算 的 必要 条 件 。 

已 有 试验 研究 表明 冻 土 中 的 未 冻 水 主要 受 土 
质 、 外 界 条 件 ( 包 括 温度 和 压力 ) .盐分 以 及 冻 融 过 
程 等 因素 的 影响 ””: 。 水 分 冻结 随 着 土 体 颗粒 粒 径 
减 小 含 盐 量 增 加 液体 饱和 度 减 小 . 土 颗 粒 活化 比 
表面 积 增 大 而 难度 增加 …。 随 着 测试 技术 的 不 断 发 
展 ,核磁 共振 法 .时 域 反射 法 和 频 域 反射 法 等 “" 广 
泛 运用 在 未 冻 水 含量 测试 中 ,通过 信号 值 与 温度 的 
对 应 关系 回归 不 同 温 度 下 的 未 冻 水 含量 值 ,建立 未 
冻 水 含量 -温度 关系 。Anderson 和 Tice"" 通 过 大 量 
室内 试验 ,依据 试验 数据 首次 发 现 未 冻 水 含量 随 温 
度 降低 呈现 指数 函数 减 小 ,并 给 出 了 经 验 公式 。 随 
后 ,此 种 关系 得 到 了 进一步 改进 及 发 展 ""”。Clap- 
eyron 方 程 建立 了 冰 、 水 压力 与 温度 之 间 的 联系 ,被 
广泛 运用 到 水 热 模型 计算 中 ,该 方程 曾 明 了 水 
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分 冻结 主要 取决 于 冰 水 压力 的 变化 。Shoop 和 Bi- 
gl “通过 土 水 特征 曲线 获取 了 冻 融 过 程 中 的 水 分 
冻结 曲线 并 预测 了 未 冻 水 含量 随 温度 变化 规律 。 
Zhang 等 ”考虑 了 土 孔 院 中 冰 对 水 的 吸附 作用 , 改 
进 了 一 般 的 土 体 冻结 曲线 。 基 于 毛细 管 压力 理论 ， 
Dall’ Amico 等 "通过 将 Clapeyron 方程 和 源 于 van 
Genuchten 模型 的 土 -水 特征 曲线 相 结合 ,建立 了 水 
分 冻结 曲线 与 土 水 势 的 关系 ,随后 就 有 学 者 证 明了 
此 方法 的 合理 性 ”。Sheshukov 和 Niber ”结合 简化 
后 的 Clapeyron 方 程 ,拓展 了 未 冻 水 含量 预测 公式 ， 
在 土 水 特征 曲线 中 引入 了 残余 含水 量 。 也 有 学 者 
根据 孔 际 大 小 分 布 函数 ,基于 孔 际 水 从 大 和 孔 际 到 小 
孔隙 的 冻结 顺序 ,建立 了 理想 成 冰 条 件 下 的 土 中 水 
冻结 模型 2。Hitchcock 55 ^ 基于 均匀 成 核 理论 、 
和 孔 际 堵塞 及 单一 孔 际 理论 曾 释 了 液体 在 介 孔 材料 
冻 融 浪 后 机 理 。 

多 和 孔 介 质 中 冰晶 预 融 现象 为 未 冻 水 含量 的 求 
解 提 供 了 新 的 思路 。Cash 等 “基于 预 融 理 论 建立 
了 负 温 下 等 粒 径 球 状 颗 粒 堆 积 体系 的 液态 水 比例 
组 成 ,并 将 未 冻 水 含量 简化 为 过 冷 度 的 函数 。Wet- 
aufer ^ 揭示 了 冰 - 液 接触 面 以 及 冰晶 表面 内 在 联 
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系 , 预 融 行为 明显 受 介质 中 杂质 的 影响 。Hendrik 和 yee MEM 


Wettlaufer“ 考 虑 了 多 孔 介 质 中 杂质 浓度 以 及 固体 
颗粒 表面 电荷 密度 的 影响 ,建立 了 理想 等 颗粒 体系 
预 融 理论 。 然 而 ,由 于 土 中 土 颗粒 粒 径 范 围 分 布 较 
广泛 , 预 融 理论 在 土 中 并 未 得 到 广泛 的 应 用 。 本 文 
在 已 有 多 孔 介 质 预 融 理 论 的 基础 上 ,提出 针对 土 体 
的 体积 未 冻 水 含量 计算 方法 ,并 分 析 了 杂质 浓度 以 
及 颗粒 粒 径 大 小 对 未 冻 水 含量 的 影响 。 


1 预 融 理论 


在 一 个 理想 等 半径 球状 固体 堆积 体系 中 ,0 °C 
以 上 纯 水 充满 固体 颗粒 间隙 ,假定 固体 颗粒 为 均 质 
的 、 各 向 同 性 且 表 面 光滑 , 且 在 超低温 下 孔 际 水 全 
部 相 变 为 冰晶 。 发 生 预 融 ( 冰 曲 开 始 融 化 ) 之 前 ,所 
有 间隙 中 充满 冰 品 , 当 温 度 上 升 后 , 冰 品 开始 融化 ， 
水 分 的 存在 形式 有 薄膜 水 和 间 际 水 ,如 图 1 所 示 。 
注 腊 水 主要 存在 土 颗粒 表面 和 冰 品 接触 面 ,而 间 际 
水 主要 存在 于 冰 蝇 与 土 颗粒 所 包围 的 间隙 中 ,其 存 
在 温度 由 冰 - 液 弯曲 界面 决定 。 


国 冰晶 体 mE DUK EE 水 膜 Emm) 土 颗粒 
图 1 多 孔 介质 中 液态 水 存在 形式 


Fig. 1 Existence form of liquid water in porous medium 


基于 多 孔 介 质 预 融 理论 ,Cash 等 “首次 建立 了 
两 种 和 常见 堆积 ,简单 立方 堆积 (Simple Cubic Pack- 
ing,SC) 和 立方 最 密 堆 积 (Face-Centered Cubic Pack- 
ing, FCC) ,其 颗粒 排列 方式 如 图 2 所 示 。 堆 积 体系 
中 液态 水 的 比例 与 薄膜 厚度 、 颗 粒 半 径 以 及 冰 - 液 
接触 面 半 径 之 间 的 关系 ,如 公式 (1) 所 示 。 


3 34-7) 1 3_r 9_r 
f(80)- 5, f, (SOM, + 5 pis? GT S r 
(la) 


f(FCO = 3 fFCOd, + 42(48 - 7/2) d, + 


1b 

342, P Dor EM 
T + 

2 R 2R? 


图 2 等 粒 径 颗 粒 排 列 方 式 


Fig.2 Arrangement of uniform size particle 


式 中 :f(SC) 和 (FCC) 分 别 表 示 SC 和 FCC 排列 下 
颗粒 占 总 体积 的 百分比 (填充 率 );f,(SC)=0.5236™; 
/FCC)=0.7405"。R 表 示 颗 粒 半径 ;d, 表 示 冰 颗粒 
和 固体 颗粒 之 间 的 水 膜 厚度 ;ds 表示 冰 颗 粒 之 间 的 
水 膜 厚度 ;所 是 球体 颗粒 占 总 体积 的 比例 。r 表 示 冰 
液 接触 面 的 半径 ;其 大 小 会 影响 冻结 温度 ;可 通过 
Gibbs-Thomson 关系 确定 1。 


若 不 考虑 系统 杂质 及 表面 带电 时 ,国体 颗粒 表 
Tai ERE BE d n] HAS (2 TCI, 
AQ, Y 

du) ii 


式 中 :人 7 表示 过 冷 度 , 即 初始 冻结 温度 降低 值 (与 
273.15 开 的 偏差 值 )。L 表 示 冰 水 相 变 潜 热 (J ,mol )。 
Pp! 是 孔 际 溶液 摩尔 密度 (mol:m”)。L 是 纯 水 冰点 温 
度 ,273.15 Ko Aud Hamaker 常数, 不 同 物质 接触 
形式 会 影响 Hamaker 常数 的 取 值 , 故 dj, 和 dd, 存 在 一 
定 差异 。 

对 于 土 而 言 ,系统 中 含有 杂质 (盐分 ) ,况且 很 
多 类 型 的 土 颗粒 表面 带 有 电荷 , 故 多 孔 介 质 体系 中 
应 该 考虑 杂质 和 颗粒 带电 量 的 影响 。Hendrik 和 
Wettlaufer ^ 推导 了 含 人 杂质 .带电 堆积 等 粒 径 颗粒 系 
统 中 的 水 膜 厚 度 与 过 冷 度 之 间 的 关系 ,如 公式 (3) 
所 示 : 


La ET ls EE 
—ATZRT, 一 zt 1- e 
T, pid DOTd E&P; cJn d 


(3) 
式 中 :ce 是 常数 ,可 通过 公式 (4) 计 算 。 
_ | eN, 
c= foe (4) 


式 中 :mm 是 土 颗粒 表面 的 杂质 浓度 ;so 是 真空 介 电 常 
数 ;s 是 液态 水 的 相对 介 电 常数 ;N\ 是 阿 伏 伽 德 罗 和 党 
BL; T,2273.15 Kie 是 电子 电荷 站 是 玻 尔 效 曼 常数 ;9， 
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是 固定 颗粒 表面 的 电荷 密度 , 孔 辽 溶液 摩尔 密度 po， 
在 计算 中 可 认为 是 常数 ;R 是 气体 常数 。 公 式 (3) 厂 
边 第 一 项 表示 杂质 浓度 影响 ,第 二 项 表示 色散 力 影 
响 下 的 界面 自由 能 ,第 三 项 表示 颗粒 表面 电荷 相关 
的 界面 自由 能 。 

以 上 理论 表明 等 颗粒 多 孔 介 质 中 水 分 冻结 问 
题 得 到 了 较 完 美的 解决 ,然而 对 于 不 同 粒 径 组 成 的 
土 中 水 分 冻结 依然 缺乏 有 效 的 方法 。 


2 研究 方法 


由 于 土 体 颗粒 分 布 范围 较为 广泛 ,例如 粉 质 黏 
土 通常 从 几 个 微米 到 几 个 毫米 , 故 现 有 预 融 理 论 在 
土 中 应 用 受到 限制 。 假 定 土 颗粒 为 球状 颗粒 ,忽略 
表面 粗糙 度 , 视 表面 为 理想 光滑 球面 。 基 于 土 体 颗 
粒 级 配 曲线 , 取 足 够 小 的 颗粒 粒 径 范 于 ,通常 试验 
测定 的 也 是 某 一 粒 径 范 围 尺 土 颗粒 质量 (体积 ) 占 总 
颗粒 质量 (体积 ) 比 例 w,, 颗 粒 分 布 及 粒 径 划 分 如 图 3 
所 示 。 根 据 参 考 文献 [28-29] ,黄土 和 砂 士 的 颗粒 
分 布 曲 线 如 图 4 所 示 。 


体积 分 数 /mm 


" " 1 " 
R; 10 100 1000 
颗粒 粒 径 /um 


图 3 粉 质 黏土 粒 径 分 布 
Fig. 3 Soilparticle distribution of silty clay 


在 第 ;个 粒 径 范围 内 ,认为 粒 径 均匀 分 布 , 土 颗 
粒 半径 均 为 尺 , 若 不 考虑 土 的 结构 性 差异 ,理想 认为 
土 体 是 以 上 两 种 颗粒 堆积 , 则 可 利用 联合 公式 (1) 
和 公式 (3) 求 解 该 粒 径 下 土 中 未 冻 水 含量 。 则 总 的 
未 冻 水 含量 可 用 公式 (5) 求 解 。 


f8O= Yo, 


3 34-m) 1 ae oe 
E f,80d,* Gedy I R RT z 
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0.14 上 — Pt 
012 上 上 M 
g 0.10 F Vul 
t V 
3 0.08 F 
S 0.06 + 
£ 
0.04 F 
0.02 F 
0.00 r= — ia aa 
0.1 1 10 100 1000 
颗粒 粒 径 /hm 
图 4 黄土 和 砂 土 粒 径 分 布 
Fig. 4 Soil particle distribution of loess and sand 
f(FCO)- > a, 
i=l 
3 1 342 r? 42r? 
= d, * 42 -mT/2 * 十 
gf FCOR, + 28 m) gd, y mrs ey 
(5b) 


式 中 :x 表示 土 颗粒 半径 为 RR 情况 下 的 冰 液 接触 面 
的 半径 。 然 而 ,公式 (5) 计 算 较 为 繁琐 ,利用 了 相同 
方法 重复 求解 不 同 粒 径 下 的 未 冻 水 含量 。 针 对 此 
问题 ,本 文 根 据 不 同 粒 径 土 颗 粒 的 体积 占 比 ,提出 
等 效 粒 径 的 方法 ,具体 计算 如 公式 (6) 所 示 。 

R= Yo, (6) 


式 中 :及 表示 土 颗粒 等 效 粒 径 。 

土 中 土 颗粒 排列 处 于 SC 和 FCC 之 间 , 为 了 进 一 
步 简 化 计算 ,本 文采 用 平均 堆积 颗粒 排列 , 即 取 胡 
(SC) 和 名 (FCC) 的 平均 值 。 结 合 公 式 (5) 和 公式 
(6) , 土 中 液态 含水 量 可 用 公式 (7) 计 算 。 

f=1.896 40.275 +7.38 — (m 
R, R, R? 


其 中 冰 - 液 接触 面 弯 曲 半 径 r 可 通过 公式 (8) 
ipie. 


pe Ya (8) 


La Nin 
p{ gear RT z] 


XP ip AE UK ra B EE OB FE 5 yo Ae VCI M A HH Be o 
对 于 冰 水 相 变 潜 热 ,在 温度 较 小 范围 内 变化 时 , 忽 
略 其 变化 , 取 厂 = 6000 J+ mol". 

随 着 温度 升 高 无 限 接近 冰点 时 ,假定 冰晶 处 于 
孔 辽 中心 处 ,无 限 小 可 近似 认为 一 点 ,如 网 5 所 示 ， 
此 时 冰晶 具有 最 大 水 膜 厚度 。 则 可 依据 颗粒 排列 
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E 冰晶 体 Em 孔隙 水 国土 颗粒 
(a) SC 模型 最 大 水 膜 厚 度 (b) FCC 模 型 最 大 水 膜 厚度 


图 5 最 大 水 膜 厚度 示意 图 


Fig.5 Schematic diagram of maximum thickness of water film 
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几何 关系 建立 不 同 土 颗粒 堆积 下 的 孔隙 初始 体积 
度 与 水 膜 表 面 杂质 浓度 之 间 的 关联 ,如 公式 (9) 所 示 。 
n; (SC) 


eC) = — > — 9a 

69 N(42 - DR, 
2n, (FCO) 

o ECO: — 2 9b 

eoo N(43 - DR, i 


式 中 :NN 为 单位 分 子 杂 质 电解 离子 数 , 对 氧化 钠 而 
言 ,N=2。 

为 了 方便 计算 土 中 初始 浓度 ,把 孔 际 水 的 初始 
浓度 co 考虑 成 SC 和 FCC 模 型 的 平均 浓度 , 则 杂质 浓 
度 可 用 公式 (10) 计 算 : 


_ ns (F CC) + n, (8C) às 
2 


假定 初始 浓度 co(SC)=coFCC)=co, 将 公式 (9) 变 
形 后 带 入 公式 (10) ,可 得 到 土 中 平均 杂质 浓度 计算 
公式 , 见 公式 (11)。 
n,, ~ 0.ANR,c, (11) 
对 于 固体 颗粒 ,Hamaker 常 数 通常 取 9.86x1072 
J ;对 于 冰晶 接触 面 而 言 ,表面 电荷 密度 可 取 0.01 
Cn 以 及 Hamaker 常 数 为 -3.3x10-2 J?" ; R=8.314 
J: (mol - K)'?;2,28.85419x10? F* m^; 6=80.4 F- 
m? ; e=1.602x10® C2”); k=1.308x10 J-K"??; N= 
6.002 x 10” mol "7 ; p=5.555 x 10* mol «ni ^"? ; p= 
5.094x10* mol: m? ?* ;y,20.03 J- m? ?* , 
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31 土 中 水 膜 厚度 与 未 冻 水 含量 变化 
一 般 土 颗粒 表面 带 有 电荷 ,以 青藏 高 原 粉 质 忒 


土 为 研究 对 象 , 取 0.1 Com?) ,该 类 土 的 初始 冻结 
温度 为 -0.2 C^", 已 有 试验 研究 表明 该 冰点 温度 归 
答 于 土 中 离子 浓度 ,基于 溶液 冰点 计算 公式 3 可 
反 算 出 初始 离子 浓度 约 0.07 ml L's 假定 土 中 盐 
分 类 型 为 典型 氧化 钠 盐 ,N=2 ,根据 公式 (9a) 估 算 n 
J 4.5x10? pmol*:m”。 以 图 3 粉 质 茜 土 颗粒 级 配 曲 线 
为 例 , 土 颗粒 等 效 半径 为 8 jym。 依 据 公 式 (3) 可 建 
立 土 颗粒 水 膜 厚度 与 冻结 温度 降低 之 间 的 关系 , 计 
算 结果 如 图 6 所 示 。 图 中 和 斜 线 表示 过 冷 度 随 水 膜 厚 
度 的 变化 规律 曲线 , 线 上 点 五 Ko RR AP IRIURE RE 
时 过 冷 度 值 。 

从 图 6 可 以 看 出 , 土 中 水 膜 厚 度 与 温度 呈现 出 
线性 关系 ,要 使 水 膜 厚度 变化 小 , 则 要 更 大 的 过 冷 
度 或 者 更 低 的 温度 。K, 和 两 个 点 表明 要 使 水 膜 
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图 6 土 中 过 冷 度 与 水 膜 厚度 关系 
Fig. 6 Curve of relationship between degree of 


supercooling and thickness of water film in soils 
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厚度 变 成 22 nm. 15 nm, 则 需要 过 冷 度 分 别 达到 20 
或 30K。 知 土 的 初始 冻结 温度 为 0% , 土 体温 度 
从 -20 *CZE 4L 28 -30 % ,水 膜 厚度 仅 变化 7nm。 该 规 
律 表 明 在 -30 WC 以 下 土 中 水 膜 厚度 变化 很 小 ,一般 
冻 土 地 区 土 中 温度 很 难 使 水 膜 厚度 变 为 几 个 纳 
米 。 故 计算 土 中 未 冻结 含水 量 ,通常 考虑 -30 CL 
内 变化 。 

类 似 地 ,冰晶 间 的 水 膜 厚度 与 过 冷 度 可 依据 公 
式 (3) 建 立 相 互 关联 ,弯曲 冰 - 液 界面 的 半径 可 通过 
公式 (8) 求 解 。 将 水 膜 厚 度 与 弯曲 液 面 半 径 的 值 代 
入 公式 (7) 即 可 求解 不 同 过 冷 度 下 的 土 中 未 冻 水 含 
量 , 当 取 等 效 半径 的 0.28 倍 计算 时 ,计算 值 与 吻合 
数据 较 好 ,计算 结果 如 图 7 所 示 。 图 中 试验 数据 均 
来 自 青 藏 高 原 粉 质 黏 士 "5235 53 ,可 以 发 现 试验 值 与 
计算 值 吻合 较 好 ,尤其 在 剧烈 冰 水 相 变 阶 段 (-2 ~ 
0 C)。 图 7 表明 改进 的 公式 可 以 求解 负 温 下 粉 质 夭 
土 中 未 冻 水 含量 。 


10r 
o 试验 数据 Pa 
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图 7 粉 质 黏土 未 冻 水 含量 计算 值 与 试验 值 对 比 


Fig. 7 Comparison between the calculated and experimental 


values of unfrozen water content in silty clay 
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为 了 进一步 验证 模型 的 适用 性 ,选取 不 同类 型 
黄土 .人 砂 土 及 条 土 试验 值 进行 对 比 。 不 同类 型 土 的 
颗粒 表面 电 葆 量 差异 较 小 ", 故 在 计算 中 均 取 0.1 
Cm”。 通 过 公式 (6) 可 计算 出 黄土 和 砂 土 的 等 效 粒 
径 分 别 为 10 jum 和 56 km。 黏土 的 等 效 粒 径 大 约 为 
4 pm?! 

不 考虑 土 中 水 过 冷 度 ,表面 杂质 浓度 与 之 前 粉 
Soa hi “LF, ,同样 对 等 效 半 径 进 行 缩小 , 黏 士 BE 
土 . 砂 土 分 别 以 等 效 半径 的 0.24、0.3、0.36 倍 计算 ， 
未 冻 水 试验 数据 均 来 源 于 参考 文献 “* 绰 , 计 算 结 
果 如 图 8 所 示 。 从 图 8a 可 以 看 出 ,计算 值 能 较 好 反 
映 未 冻 水 变化 规律 , 当 温 度 较 低 时 , 即 A7> 10 K, fit 
测 值 比试 验 值 较 高 。 对 于 砂 土 而 言 , 计 算 值 与 试验 
值 吻合 较 好 ,未 冻 水 体积 含量 迅速 减 小 并 趋 于 稳 
定 , 当 A7T>5K 时 ,未 冻 水 比例 基本 不 变 。 而 对 于 条 
土 ,虽然 计算 值 反 映 了 变化 趋势 ,但 误差 较 大 ,尤其 
在 开始 冻结 温度 段 ,其 主要 原因 是 条 土 颗粒 的 形状 
星 片 状 , 与 假定 的 球形 差距 较 大 。 

3.2 水 膜 厚度 影响 因素 

不 同 物理 化 学 界面 能 对 水 膜 厚度 的 影响 如 图 9 
所 示 。 色 散 力 其 与 Hamaker 常数 相关 ,对 水 膜 的 影 
响 与 杂质 浓度 没有 关系 , 随 着 水 膜 厚度 减 小 ,过 冷 
度 呈 线性 增 大 , 当 水 膜 减 小 至 几 个 纳米 级 时 ,其 逐 
渐 发 挥 作 用 。 杂 质 浓 度 对 水 膜 厚 度 的 变化 影响 很 
大 , 随 着 浓度 增加 ,相同 过 冷 度 下 水 膜 厚度 减 小 越 
来 越 困 难 ,说 明 水 分 冻结 难度 增 大 。 土 颗粒 表面 电 
答 密 度 对 水 膜 厚 度 的 影响 极 易 受 杂 质 浓度 的 影响 ， 
随 着 浓度 增 大 , 土 颗粒 表面 电 葵 对 水 膜 的 影响 越 来 
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图 8 不 同 土 类 型 计算 值 与 试验 值 对 比 


Fig. 8 Comparison between the calculated and experimental values of different soils 
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图 9 不 同 杂 质 浓度 及 物理 化 学 界面 能 对 水 膜 厚度 影响 
Fig.9 Effect of different impurity concentration and physical 


chemical interface energy 


越 小 ;相反 , HAR CT BE Da) 0.01 umol: m?, EI 
To TT HEL qup FRE HT EFL B PK BERE E BS ME HEC 
杂质 浓度 本 身 的 影响 。 当 杂质 浓度 大 于 1 pmol*m?, 
表面 电荷 密度 对 水 膜 的 影响 非常 小 ,在 计算 中 可 以 
不 考虑 。 此 外 ,依据 公式 (3) , 当 颗 粒 表面 电荷 密度 
减弱 时 ,其 对 水 膜 的 影响 也 随 之 减 小 。 
3.3 薄膜 水 和 间隙 水 变化 规律 及 影响 因素 

对 于 粉 质 黏 士 ,颗粒 表面 水 膜 ,冰晶 间 水 膜 以 
及 间 际 水 所 引起 未 冻 水 含量 的 变化 如 图 10 所 示 。 
可 以 看 出 , 土 中 水 分 冻结 主要 由 土 颗粒 表面 水 膜 厚 
度 变 化 来 决定 。 冰 唱 间 水 膜 的 影响 相 比 土 颗粒 表 
面 水 膜 而 言 较 小 ,两 者 计算 值 之 间 差 了 一 个 数量 
级 。 而 土 中 间隙 水 对 未 冻 水 含量 的 影响 最 小 ,通过 
公式 (8) 得 出 , 较 高 的 初始 浓度 会 减 小 冰 - 液 弯曲 界 
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105 D 基于 土 颗粒 表 H 


图 10 水 膜 及 间隙 水 对 未 冻 水 含量 计算 值 的 贡献 


Fig. 10 Construction of water film and interstitial water to 


calculation of unfrozen water content 


面 的 半径 ,导致 7 及 .的 平方 减 小 很 多 ,致使 液态 含水 
量 很 小 。 

当 土 中 杂质 浓度 变化 或 土 颗粒 粒 径 变化 ,诸如 
黏土 比 粉 质 黏 士 的 颗粒 要 细 ,在 此 假定 不 同 浓度 及 
粒 径 条 件 ,未 冻 水 含量 计算 结果 如 图 11 所 示 。 可 以 
发 现在 浓度 较 高 的 条 件 下 (图 1la) , 随 着 士 的 等 效 
颗粒 粒 径 减 小 ,未 冻 水 含量 增 大 ,表明 水 在 细 颗 粒 
土 中 更 难 冻结 。 在 浓度 较 低 情况 下 (图 11b), 未 冻 
水 含量 变化 表现 出 相似 的 规律 ,但 基于 间隙 水 计算 
得 到 的 未 冻 水 含量 所 发 挥 的 作用 随 着 颗粒 减 小 愈 
发 明显 ,在 很 大 温度 范围 内 ,超过 了 基于 冰晶 间 水 
膜 厚 度 计算 的 未 冻结 含水 量 值 。 此 外 ,计算 结果 表 
明 土 颗粒 等 效 粒 径 越 细 , 初 始 冻 结 温度 较 低 , 低 浓 
度 细 颗粒 的 冰 - 液 弯曲 界面 在 水 分 初始 冻结 中 占 主 


wp Of — 基于 间 阶 水 弯曲 冰 - 液 界面 计算 
E 一 基于 冰晶 界面 水 膜 厚度 计算 
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图 11 浓度 及 颗粒 半径 对 未 冻 水 含量 的 影响 


Fig. 11 Effect of concentration and particle radius on unfrozen water content 
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导 地 位 。 易 受 杂质 浓度 的 影响 , 随 着 浓度 增 大 , 土 颗粒 表面 


4 讨论 


经 典 的 指数 或 星 函 数 关 系 * 在 水 分 剧烈 相 变 阶 
段 ( 初 始 冻 结 温 度 附近 ) 存 在 很 大 的 误差 ,因为 指数 
函数 在 冰点 0 不 收敛 , 且 模 型 中 参数 随 不 同 初始 
含水 量变 化 较 大 。 结 合 土 水 特征 曲线 和 冻结 特征 
曲线 的 预测 模型 ”很 大 程度 上 依赖 于 初始 含水 量 
和 残余 含水 量 ,上 且 模 型 中 涉及 的 参数 较 多 。 这 些 未 
冻 水 含量 模型 中 的 参数 没有 明确 的 物理 意义 ,而 本 
文 改进 的 模型 物理 参数 明确 , 且 考 虑 了 土 体 颗粒 大 
小 .杂质 浓度 、 颗 粒 表面 带电 特性 ,并 涉及 了 水 分 冻 
结 的 物理 化 学 微观 机 理 。 

Dash 等 ”通过 预 融 理 论 计算 了 等 粒 径 聚 茶 乙 
烯 粉 末 及 石墨 化 炭 黑 粉 中 未 冻 水 含量 的 分 布 ,发 现 
当 半 径 分 别 缩小 0.15 和 0.4 时 计算 效果 与 试验 值 吻 
合 较 好 。 土 颗粒 实际 排列 非常 复杂 , 且 颗 粒 形状 并 
非 全 是 球状 体 , 且 士 颗 粒 表面 呈现 出 一 定 粗糙 度 。 
本 文 在 计算 过 程 中 采用 了 等 效 粒 径 , 当 粉 质 黏 士 、 
黄土 . 砂 士 分 别 以 等 效 半径 的 0.28 .0.3.0.36 倍 计算 
时 ,计算 值 与 试验 值 吻合 较 好 。 而 对 于 细 颗 粒 土 黏 
土 ,计算 值 与 试验 值 差别 较 大 。 此 外 ,在 模型 中 未 
考虑 细 颗 粒 土 中 特殊 结构 ,大 部 分 黏土 中 存在 如 蜂 
窝 结 构 .名状 结构 等 特殊 形式 ,与 自然 排列 下 的 土 
颗粒 分 布 有 一 定 差距 ,会 影响 模型 计算 精度 。 同 
时 ,计算 中 缩小 半径 的 处 理 方 式 还 未 有 详细 阐述 。 
所 以 颗粒 形状 .颗粒 表面 粗糙 度 . 土 颗 粒 的 实际 排 
列 和 土 的 结构 类 型 等 会 作为 水 分 微观 冻结 方面 后 
期 研究 内 容 , 进 而 开展 进一步 研究 。 


5 结论 


本 文 基于 多 孔 介 质 预 融 理论 ,开展 了 土 中 未 冻 
水 体积 含量 随 温度 变化 研究 ,提出 了 其 计算 方法 ， 
并 分 析 讨 论 了 杂质 浓度 和 粒 径 大 小 对 水 膜 厚度 及 
未 冻 水 含量 的 影响 ,具体 结论 如 下 : 

(1) 对 于 较 大 颗粒 土 ,未 冻 水 含量 模型 预测 效 
RR. HWER E 黄土、 砂 士 等 效 粒 径 分 别 缩 
小 0.28 .0.3 .0.36 倍 时 ,未 冻 水 含量 计算 结果 与 试验 
值 接近 ,模型 可 较 好 地 预测 土 中 不 同 过 冷 度 下 液态 
水 含量 。 而 对 于 黏土 ,未 冻 水 含量 预测 误差 较 大 。 

(2) 土 颗 粒 表面 电荷 密度 对 水 膜 厚 度 的 影响 极 


电荷 对 水 膜 的 影响 越 来 越 小 。 随 着 杂质 浓度 增加 ， 
相同 过 冷 度 下 水 膜 厚度 增加 ,水 分 冻结 难度 增 大 。 

(3) 土 中 未 冻 水 含量 计算 可 以 不 考虑 色散 力 的 
影响 ,水 分 冻结 主要 由 土 颗粒 表面 水 膜 厚度 变化 来 
决定 。 当 等 效 颗粒 减 小 时 ,通过 间 际 水 求解 的 未 冻 
水 体积 比 逐 渐 增 大 ,并且 在 较 低 浓度 时 发 挥 更 大 
作用 。 

(4) 当 土 颗粒 表面 杂质 浓度 小 于 1 pmol s m^, 
土 中 水 分 冻结 温度 由 冰 - 液 接触 面 半 径 决 定 , 该 半 
径 随 等 效 粒 径 减 小 而 减 小 , 且 半 径 越 小 ,冻结 温度 
越 低 。 
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Variation of water film thickness in soil and prediction 
method of unfrozen water content 


WAN Xusheng, YAN Mengyu, LU Jianguo， YAN Zhongrui 


(School of Civil Engineering and Geomaties, Southwest Petroleum University, Chengdu 610500, Sichuan, China) 


Abstract: The variation of unfrozen water content directly influences the thermodynamics and deformation 
properties of frozen soils, and it is also the key condition for water-heat coupled simulations. In order to study 
this variation, the distribution of soil particle size was used to an equivalent particle size, and soil was simplified 
to the equivalent particle size spherical packing system. Based on the premelting theory in porous medium, the 
calculation method of unfrozen water content in soil was put forward with the considered average packing system 
of simple cubic packing and cubic close packing, and the accuracy was verified by the experimental data. 
Moreover, the effect of impurity density and the equivalent particle size on water film thickness and unfrozen 
water content was analyzed. The results show that the thickness of water film calculated using the surface charge 
density of soil particles was easily affected by impurity concentration, and the effect of surface charge density on 
water film was progressively reduced as the impurity concentration increased. The liquid water fraction was 
mainly determined by the variation of water film thickness on the surface of soil particles, and the contribution of 
interstitial water to the total liquid fraction increased as the equivalent particle size decreased, especially at low 
impurity concentrations. Moreover, the model of unfrozen water content produced better predictions in larger 
particle soils. The calculated values of volumetric unfrozen water were close to the experimental values when the 
equivalent particle sizes of silty clay, loess, and sand were reduced by 0.28, 0.3, and 0.36, respectively. 


Keywords: premelting; unfrozen water content; equivalent particle size; impurity density; surface charge density 


